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Abstract 

We estimate the flux of the GRB neutrino background and compute the 
event rate at SK and TITAND in the collapsar model, assuming that GRB 
formation rate is proportional to the star formation rate. We find that the flux 
and the event rate depend sensitively on the mass-accretion rate. Although 
the detection of signals from GRBs seems to be difficult by SK, we find that 
we can detect them by TITAND (optimistically ~ 20 events per year in the 
energy range En e > 50MeV). Therefore, we show that we will obtain the 
informations on the mass-accretion rate of collapsar by observing the GRB 
neutrino background, which, in turn, should give us informations on the total 
explosion energy of GRBs. 
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I. INTRODUCTION 



Gamma-Ray Burst (GRB) is one of the most mysterious phenomena in the universe. 
Although there was a great revolution in our understanding on GRBs in 1997 that some 
of them were proved to be of extragalactic origins QXJ] , there are still a lot of questions and 
puzzles left on the phenomena. That is why GRBs have been the fascinating targets for 
astronomers. In fact, there are a lot of interesting and valuable discussions under debate. 
For example, we have not known what is the origin of GRB. Although this problem has 
never been settled and still controversial, there are some reports that suggest the physical 
connections between GRBs and supernovae (SNe) [Q. If the SN/GRB connection is correct, 
GRBs will be born from the death of massive stars. Then the GRB event rate should trace 
the star formation history in the universe 0. 

The central engine of a GRB is an important but difficult problem to be solved. Although 
a fireball model || is a promising and attractive one, we do not know how the initial condition 
is. We merely gave it artificially. In the fireball model, since the ratio of baryon number 
to photon number is required to be extremely small as an initial condition, only a peculiar 
environment could meet such a condition. As long as we believe the GRB/SN connection, 
however, a collapsar (failed supernova) model || will be one of the convincing candidates 
as the central engine of a GRB. In the model, much energy is generated by the neutrino 
heating from an accretion disk surrounding a black hole which is made by the core collapse 
of a massive star. 

On the other hand, total explosion energy of a GRB is also one of the most important 
problems in the physics of GRBs. Observationally it is well known that there is a large 
dispersion of total energy of gamma-rays emitted from GRBs ||. Although the dispersion 
may become small when we consider the beaming effects 0, after all we do not know 
the total explosion energy of the system that generates a GRB only by the gamma-ray 
observations.^ In the theoretical side, the total explosion energy in the collapsar model is 
also not estimated clearly, and only a rough estimate is reported as a result of the numerical 
computations ||. Therefore, to know the total explosion energy, we require the multi- 
dimension numerical hydrodynamical computations including the microphysics of neutrino 
heating in future, which is, of course, a laborious task to be performed. 

Are there observations which give us new clues to know how a GRB is produced? In 
this study, we point out that the GRB neutrino background may be detected in a large 
water Cherenkov neutrino detector such as Super-Kamiokande (SK) and TITAND which is 
proposed as a next-generation plan of the multi- megaton water Cherenkov detector |J . We 
also stress that this signal should give us valuable informations on GRBs, such as the total 
explosion energy and/or event rate of GRBs. The reason is as follows. In this study, we 
assume that the origins of GRBs are massive stars. Thus, the event rate of GRBs should 
trace the star formation history in the universe. Since we adopt the collapsar model here, we 
can assume that the central engine of a GRB is the neutrino heating from an accretion disk 
surrounding the central black hole. We show that the total number of neutrinos emitted 



1 Here, the total explosion energy means the summation of energies of photons, leptons and 
baryons emitted from an origin of a GRB except for neutrinos. 
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from a GRB is very sensitive to the mass-accretion rate M which is difficult to estimate 
only from the existing theoretical models. However, we can obtain the informations of 
the mass-accretion rate by observing the GRB neutrino background. We also derive an 
approximate relationship between the total explosion energy and the mass-accretion rate of 
a GRB. Namely, although we cannot determine the total explosion energy only from the 
gamma-ray observations, we can obtain informations of it by observing the GRB neutrino 
background. That is an important clue to understand GRBs. In this study we find that 
the GRB neutrino background may be detected within ~ 10 yrs by TITAND as long as the 
average mass-accretion rate among GRBs is M > a few M® s -1 , and the probability that 
one collapsar generates a GRB is / ~ 0.5 - 1.0. 

Here we outline the organization of this paper. In section H, we explain the formulations 
in this study. Results are shown in section |T|. Discussions are given in section [TV]. Summary 
and conclusion are presented in section [V]. 



II. FORMULATION 

A. GRB formation history 

First we assume that the rate of GRBs traces the global star formation history of the 
universe, -Rgrb(- 2 ) oc -Rsfr- Although a number of workers have modeled the expected 
evolution of the cosmic SFR with redshift, there are some uncertainties (in particular, at 
high redshift z > 1), and further observations are required to determine the cosmic SFR. In 
this study, we use the three different parameterizations of the global star formation rate per 
comoving volume in an Einstein-de Sitter universe to take into account the uncertainty of 
the cosmic SFR which is derived by Ref. ||. They are 

R SF1 (z) = 0.3/i 65 M^Mpc" 3 , (1) 

exp(3.8z) + 45 

R SF2 (z) = 0.15^5 ^i 3 ; 4 l } 99 M^r^Mpc" 3 , (2) 
exp(3.4;z) + 22 

and 

^ , n , exp(3.05z — 0.4) , 

fe < 2 > = °- 2k '» exp(2.93,) + 15 M ^ M ^ < 3 > 

where h e5 = H /65 km s -1 Mpc -1 . 

Assuming that massive stars with masses larger than ~ 35M Q explode as GRBs |i~0" 



GRB-rate density -Rgrb can then be estimated by multiplying the selected SFR by the 
coefficient 




x , g ^^ = 1.5xl0-M^)M-i, (4) 



Jg dm 0(m) 
Jq 1 ^ 5 dm m<fi(m) 

where 4>(m) is the initial mass function (IMF) and m is the stellar mass in solar units. In this 
study, the Salpeter's IMF (<f>(m) oc mr 2 - 35 ) is adopted throughout. Here we assumed that one 
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tenth of the collapsars whose mass range is in (35 - 125)M are accompanied with GRBs (= 
/). The validity of this assumption and the detailed discussions are presented in section |V|. 
We label these models as GRBR1, GRBR2, and GRBR3, respectively. For comparison, we 
use the SN-rate density i?sN by multiplying the selected SFR by the coefficient 

£ d ™*M = 1.2 x lO-M- (5) 
Jo 4 dm m<j){m) 

As for the energy spectrum of V e of SNe, the Fermi-Dirac distribution with zero chemical 
potential is used (see Figs |9|- |11|), which have the same luminosity with the numerical 
model Jllj . 

In the flat universe, i.e., when the density parameter is Q = Qm + — 1 with the 
non-relativistic matter density (Qm ~ 0.3) and the A term (Q\ ~ 0.7), the star formation 
rate is given by |3f 



'((1 + n M z))(l + zf - n A (2z + z 2 ) 
R S f(z; Qm-^a, h 65 ) = ^65 J 7— — rT/9 Rsf(z; 1, 0, 1), (6) 

(i + zy /z 

where z is redshift. In Fig. [I], we show the GRB rates as a function of the redshift. 



B. Accretion disk model 

1. Energy Spectrum 

In the collapsar model, neutrinos are mainly generated at the inner-most region of an 
accretion disk which is formed around the central black hole. Thus we estimate the energy 
spectrum of neutrinos using an analytical shape of the accretion disk. Popham et al. (1999) 
proposed an accretion-disk model including the effects of general relativity as a collapsar 
model [P^j . The model in Ref. |12j and the numerically computed collapsar model are 
fitted well by the analytical shape of the accretion disk in Ref |]T2]JT3[| 2 . We use it in this 
study. Then, the density [g cm -3 ], the temperature [MeV] and the disk thickness [cm] are 
fitted as follows (see Ref. |15[] for details). 



p = 8.23 x10 s — - TTTxr-T T^TWT T^Tore; (7) 




T = 2-3 x — - TTTTFTZT 77^ TTTTVolT' (») 



2 On the other hand, in the merger models, e.g., for the merger of NS-NS binaries, the analyti- 
cal equilibrium solutions of the neutrino-cooled accretion disk and the properties of the neutrino 



emission are discussed in Ref. [14|. 
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/ v 0.9 / \ -0.0183 

# = 5.8xl0 6 ( — ] ( ^— r ) ™ ™, (9) 

V3M y ^o.iMos- 1 ; ^ ^ 1 + ^y 

where M is the mass of the central black hole, r is the radial coordinate, and r s (= 10 7 cm) 
is the core radius, respectively. Note that the Schwarzshild radius is ~ 8.862km(M/3M Q ). 
In Ref. |p.2|| , they assumed that the mass of the central black hole is constant. In this study, 
however, to estimate the total energy of neutrinos we consider the effects that the mass of 
the central black hole becomes larger with time by accretions. 

When we know the density and the temperature in the neutrino emitting region, we can 
estimate the emissivity of neutrinos (cm _3 sec _1 MeV _1 ) in the accretion disk. The total 
emissivity of V e is represented by two terms: 

d 2 n v d 2 n% N d 2 n^ +e ~ , , 

— = — H - — (10) 

dtdE Ve dtdE Ve dtdE Ve ' v ; 

where Ejj e is energy of V e . For n + e + — > p + V e , the emissivity of V e is represented by 



dtdEu. (EvJ = 2^ (1 + 3 &K^.V(^. - Q) - K {E Ve - Q) e(Ehe _ Q)/T + v (") 
where Gf is Fermi coupling constant, qa is the axial vector coupling constant of the nucleon 



which is normalized by the experimental value of neutron lifetime r n ~ 886.7 s flgl , n n is 
number density of neutron, Q ~ 1.29 MeV, and m e is electron mass. For e~ + e~ — > u e + T7 e , 
we obtain 

(£VJ = g (c . + cm*. —'t- r (g2 - ( .?: / r )2)3/2 rf , (12, 



dtd^ v UeJ 9n^ V AJ ^e E ^J T + l Jm e /T e e + l 



where Cy = 1/2 + 2 sin 9w, Ca = 1/2, and sin 6w — 0.231 is Weinberg angle [|T6|] , and we 
assume T 3> m e . We plot them in Fig. |2| in the case that T = 5MeV and p = 10 10 g • cm"" 3 . 
Using Eqs. (0), (||) and (Q), we obtain the integrated energy spectrum of u e , 

dn , . /" A * , r , cPr^ . . , . 

■{Eh.) = dt , d v a —£-(Eh B ), (13) 



d£x7 e ^ e Jo JT>iMeV a dtdE-y^ 

where / dV a denotes the integration in the volume of the emitting region in accretion disk 
where the temperature is high enough to produce neutrinos, i.e., T > 1 MeV, / dt denotes 
the integration of time t, and At = M a /M is the duration of the accretion (M a is the mass 
of the progenitor). Here we assume that the time evolution of the mass of the central black 
hole is given by 

M = Mi + Mt, (14) 

where Mj is the initial value of M. In this study, M% is set to be 3M . In Fig. [5] we plot the 
obtained energy spectrum of V e emitted from the accretion disk in unit energy [MeV -1 ]. 
Here we define the time-dependent luminosity of v e and V e by 
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where dn Ve /dtdE Ve is the total emissivity of z/ e with their energy i^ e , and we assume it is 
approximately equal to dnn e /dtdE Ve . In addition, we define the total energy of u^s and V^s 
emitted from a collapsar, 

/■At 

Qtot = / L v (t)dt. (16) 
•/ o 

We see from Eq. (§) that as the mass of the central black hole grows, the temperature 
becomes lower, and the flux of neutrinos decreases. In Fig. |4], we show the duration (= 
At^emit) which is defined as the e-folding time of the time- dependent luminosity in Eq. (|15|). 
This almost corresponds to the timescale in which the temperature at the innermost region 
of the accretion disk becomes lower than ~ 1 MeV. Namely this will approximately reflect 
the timescale of a GRB. 

Now we investigate the dependence of the energy spectrum of anti-electron neutrino on 
the mass-accretion rate (M). The mass-accretion rate is theoretically not known while the 
numerical simulations have been performed, because there are uncertainties on the viscosity 
of the accretion disk and angular momentum of the progenitor. In Fig. [5|, we show the 
energy spectrum for the case of M = 0.01M Q s -1 , O.1M s _1 , 1M s _1 , 5M Q s _1 , and 
1OM s _1 . The energy spectrum is represented as E^ e dn/ ' dE v& [MeV] because the cross 
section of the reaction p + V e — > n + e + is proportional to the square of the energy. It is 
clearly shown that the energy spectrum becomes hard, and intensity becomes large as the 
mass-accretion rate increases. This is because the temperature at the inner region of the 
accretion disk becomes higher. Here we have to give some comments on the application 
limit for the present accretion disk model in this study. High energy electron-positron pairs 
whose energies are greater than ~ lOOMeV produce muons and charged pions, which, in 
turn, decay into low energy leptons. Such microphysics is not included in the present model. 
The influence might become important for the high accretion rate model. In addition, as the 
density becomes large, the assumption that the neutrinos are optically thin might become 



to break down ||12|| . For example, the opacity of electron neutrino (r) at the inner most 



region of the accretion disk is estimated as |TJ 

r = J ^N A padr (17) 



1.6 x 10- 9 T^ cY p w H (18) 

1.23 



., ,3M \ L5 / M 



r. 



where Na is the Avogadro's number, a is the cross section of the neutrino off a nucleon, pio 
is the density normalized by 10 10 g cm -3 . Therefore, r may comparably become unity if M 
is approximately 10M Q s" 1 , r ~ r g , and M ~ 3M . From the above considerations, models 
with high accretion rate (M ~ 1OM s _1 ) may be the extreme cases. 

In Fig. |6], we show the average neutrino luminosity which is defined by 

L u = Qtot/At^emit- (20) 
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For comparison, in the figure we also show the result of Popham et al. (1999). The un- 
certainty of their results comes from the Kerr parameter of the central black hole (a = 
- 0.5). From the figure, we find that our simple formula reproduce their results fairly well 
in wide range of mass-accretion rate. At high mass-accretion rate (M ~ 10M Q s _1 ), our 
formula may overestimate the neutrino luminosity by a factor of three. As shown below (see 
section |T|), however, we also compute the GRB neutrino background in such models as ex- 
treme cases. In future, improvements of the accretion model will enhance the predictability 
of the flux of GRB neutrino background. 

The event rates of V e expected in water Cherenkov detectors is represented by 

dR N v . „ . dn . ^ , , , 

' a pVe (E e+ ) — (E e+ ) } (21) 



dE e+ AttD 2 pUey e ' dE v< 
where E e + = Ep e — Q is the energy of the positron which is scattered through p+V e — > n + e 



G o / 

in the detector, a pVe = + 3g\)E e +J — m 2 is the cross section of the process, and D 
is the distance from the Earth to the collapsar. N p is the number of proton in the detector, 
e.g., N p ~ 1.5 X 10 33 for SK. Thus the total event rate is estimated by the integration, 



2. Estimate for the total explosion energy 



It is reported that the energy deposition rate becomes as large as ~ 10 51 erg s 1 by 



the vV pair annihilation ||12|| . In their studies, they computed the energy deposition rate 
considering the shape of the accretion disk. Their results are shown in Fig. [?| as a shaded 
region. Here the efficiency of i/ e z7 e pair annihilations is defined as e = Q e + e -/Qtot, where 
Qe+e- is the total deposited energy. Their results are fitted well as 

/ • v 1.56 

which is also shown in Fig. [7[ In this study, we use this fitting formula to estimate the 
explosion energy of a GRB. Then, we obtain the total deposited energy by Q e + e - = eQtot- In 
Fig. §, the total energy of emitted neutrinos (Qtot) and the total deposited energy (Q e _|_ e _) 
are shown as a function of the mass-accretion rate. Since the total explosion energy of 
SN1998bw is reported to be ~ 10 52 erg M, we find that it corresponds to the mass-accretion 
rate ~ (1 - 2)M s" 1 in Fig. gf 

For simplicity, we ignored the heating process from v^p^ and v T V T pair annihilations 
in this computation because the effect is much smaller than z/ e 77 e pair annihilation. The 
reason is as follows. In the electromagnetic thermal bath, f^y ^ and u T V T pairs are produced 
from e + e~ pair annihilations only through the neutral current interaction. In addition, v^p^ 
and v T V T pairs annihilate into e + e~ pairs only through the neutral current interaction. On 
the other hand v^p,, pairs are produced through both the charged and neutral currents and 
annihilate into e + e~ pairs correspondingly through the same process. Then, compared with 
the case of u e V e , the efficiency of the heating from v^V^ and u T T7 T pairs in the electromagnetic 
thermal plasma is just within 10% of the total, even if the z/77 pair production from e + e _ 
pair annihilation dominates all the other processes of neutrino production. 
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C. GRB neutrino background 



In order to get the differential number ffux of back ground neutrinos, first we compute the 
present number density of the background neutrinos per unit neutrino energy, dn(E Vc ) / dE Ve 
in Eq. ([13]). The contribution of the neutrinos emitted in the interval of the redshift z ~ z+dz 
is given as 

dN v {Eh m ) = R GKB{z)f z dz^^-{l + z)dE Ve , (24) 

where = (1 + z)E-p e . It is noted that -Rgrb is GRB rate per comoving volume, and the 
effect of the expansion of the universe is taken into account in Eq. (|24]). The Friedmann 
equation gives the relation between t and z as follows: 

dz 



H (l + z)y/(l + n M z)(l + z) 2 - n A (2z + z 2 ). (25) 

We now obtain the differential number flux of GRB, dF u (E Ve )/dE Ve , by using the relation 

dF u (E Ve ) /d^ = cdN v (E Ve ) /dEb e : 

dF v c dn{{l + z)E Ve ) dz 

TT /_ r gkb(z) tb 7=. = = = = = =t= ( 26 ) 



dEv e H Jo ' J dE^ e ,/(l + Q M z)(l + zf - VL K {2z + z 2 



III. RESULTS 

The differential number flux of GRB neutrinos is shown in Figs. ||], |10| and [IT]. As for 
the GRB formation history, GRBR1, GRBR2, and GRBR3 are assumed, respectively. It is 
found that the differential number flux is similar with each other even if the different star 
formation history is adopted. That is because the energy of neutrinos from high-z GRBs 
is redshifted as oc 1/(1 + z), and the volume containing high-z GRBs is relatively small 
oc 1/^/(1 + Qmz)(1 + z) 2 — Q\(2z + z 2 ) as given in Eq. (PB|). Thus, we show the results 
only for GRBR3 below. On the other hand, we can see clearly that the differential number 
flux becomes larger as the mass-accretion rate increases. In particular, in high accretion 
model such as M = 5M , the intensity of the GRB neutrino background is larger than 
that of SN neutrino background in the energy region E Ve > 50MeV. As discussed below, in 
such a high energy region, the neutrino background will be obscured by the noise from the 
atmospheric neutrino background (ANB). As the event number increases, however, we will 
have statistically sufficient data, and the signals from GRBs will dominate the statistical 
error of ANB. In the next section we estimate the time required to detect the signals from 
GRB (see also Table pip . Note that at the energy which is larger than 50MeV, the signals 
from GRBs dominate the noise of electrons from the decay of invisible muons because their 
spectrum has a steep cut-off at 50MeV. In the next section, we will discuss them in detail. 

The event rates at SK are shown in Fig. [TJ| in units of yr -1 MeV -1 for various mass- 
accretion rates. For event rates, increasingly we cannot discriminate among the three models 
of the star formation history because the cross section of p+V e — > n+e + is proportional to the 
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square of the neutrino energy (oc E% ), and the relative contribution of low energy neutrinos 
to the event rate is small. Therefore, irrespective of the differences among the models of 
the star formation history, we can expect the event rate of GRB neutrino background as 
a function of mass-accretion rate. Namely the informations on the average mass-accretion 
rate will be obtained from such an observation, which, in turn, should give us informations 
on the total explosion energy of a GRB (see Fig. §). They are, of course, never obtained 
only by the gamma-ray observations. 



In Fig. [L3|, we show the dependence of the event rates on the probability / that one 
collapsar generates a GRB in the case of M — 1OM s _1 . It is clear that the signals from 
GRBs dominate those from SN for E De > 30MeV in the case of / = 1.0. However, the 
noise from the atmospheric neutrinos obscures the signals of GRBs even if we adopt the 
high value of / (=1.0). In the next section, we discuss the detectability of the signals of the 
GRB neutrino background as a significant deviation from the statistical error of ANB. 

To estimate the event rates more realistically, we should consider the effects of neutrino 
oscillation. The flavor eigen states v a (a, = e, //, and r) are related with the mass eigen states 
Vi ("i" = 1, 2, and 3) by the mixing matrix U as v a = J2i U a iVi. Then, we parameterize it 
by 



u 



1 C13C12 C13S12 S13 

-C23S12 - S23S13C12 C23C12 - S23S13S12 S 23Cl3 I , (27) 
V S23S12 - C23S13C12 -S23C12 - C23S13S12 C23C13 



with Sij = sin 9ij and Cjj = cos 9ij, where we ignored the CP phase in the lepton sector. In 
general, the probability of the vacuum oscillation from a flavor v a to a flavor vp is represented 
by 



P (u a _> Vp ) = (u p \u a (t)) = J2 \Upi?\Uai\ 2 + 2ReJ2U^U; j U* m U aj e- i ^ 5 < 1 (28) 



<5m 2 

i>j 



where L is the distance from the source to the Earth (= ct), E v is the neutrino energy, 
and the mass difference is defined by Smfj = — mf. For simplicity, here we assume 
the normal hierarchy of neutrino masses, mi < 7712 < m^. The results of recent oscillation 
experiments are as follows. The solar neutrino data favor the Large Mixing Angle solution 
(0.6 < sin 2 2# 12 < 0.98), and 2 x 10~ 5 eV 2 < 5mj 2 < 4 x 10" 4 eV 2 at 3a C.L. @. The 



atmospheric neutrino data also imply the large mixing (sin 

2^3 > 0.85), and l.lxl0" 3 eV 2 < 



(5m^ 3 (~ 5m 2 13 ) < 5 x 10 _3 eV 2 at 99 % C.L. The CHOOZ reactor experiment shows that 



sin 2 2^i3 < 0.1 at 95 % C.L p0 |. Then, the above data tells us that 5m\ 2 *C 5m\ 3 ~ 5m 



2 

13- 



Adopting the above data, we can assume that the mixing matrix is approximately sim- 
plified into 

/ V2/2 V2/2 \ 

U ~ -1/2 1/2 V2/2 . (29) 
V 1/2 -1/2 y/2/2j 

Because neutrinos propagate for long distances, the second term in Eq. fl28|) is averaged and 
becomes zero. Then, we have 



9 



P(v e -> F e ) = 1/2, (30) 

-> Fe) = 1/4, (31) 

P(F r -> F e ) = 1/4. (32) 

If we choose the inverted hierarchy, i.e., 5m\ 2 -C 5?^2i — ^ m 3i5 we only interchange the roles 
of mi and in Eq. fl27f) . Namely we only interchange the first and third columns of U. 
Even then, it is remarkable that the results in Eq. (p0[), (|31~D, and (|3^) are not changed at 
all. 



In Fig. [TJ|, we show the event rate as a function of E Ve at SK in which we consider the 
effects of the neutrino oscillation for M = 10M Q . It is noted that the flux of T7 e becomes 
lower as a result of the neutrino oscillation. That is because and V T are produced only 
through the e + e~ pair annihilation. 



IV. DISCUSSIONS 

We consider the detectability of the neutrino background from GRBs. In the energy 
range E < 50MeV, the neutrino background is obscured by the events of the electrons 



which come from the decay of invisible muons [21]]. The spectrum of the electrons is given 
by dN/dE e = {G 2 F / 12'K 3 )m 2 l E 2 (3 — 4E e /m^), where m M is the muon mass, and E e is the 
electron (positron) energy. The event rate of the electrons from the invisible muons at SK 
is about ~ 1 in units of yr -1 MeV -1 in the energy range 50 > Ep e > 20MeV [|TTJ. Thus, 
signals of GRB neutrino background are smaller than those of electrons from invisible muons 
in the range. In addition, the signals from SN neutrino background will basically dominate 
those of GRB neutrino background in the energy range. On the other hand, the spectrum of 
electrons from the invisible muons has a steep cut-off at about 50 MeV because the maximum 
energy of the emitted electrons is about half of the muon mass. Thus, as mentioned above, 
the signal from GRBs becomes larger than the noise from invisible muons in the energy 
range E T7e > 50MeV. In the case of the high mass-accretion rate such as M > 5M , the 
signals of the GRB neutrino background are larger than those of SN neutrino background 
in the energy region. However, the severest noise comes from the atmospheric neutrinos. In 
Table |, we show the event rate (yr -1 ) of ANB and its statistical error in SK/TITAND in 
the energy regions E Pe = (50 - 60) MeV and E Pe > 50MeV J22]|23|. TITAND is proposed 
as a next-generation plan of the multi- megaton water Cherenkov detector ||. The first 
phase of TITAND is planed to be with 2 Mt water inside. Thus, the event rate becomes 
approximately one hundred times larger than SK. 

For comparison, we show the event rates from the GRB neutrino background at 
SK/TITAND in Table [H] for various models. The third column of the table represents 
whether we consider the effects of the neutrino oscillation or not. 'Yes' means that the 
effects of the neutrino oscillation are taken into account. It is apparent that the signals from 
GRB are smaller than the noise from ANB. As the event number increases, however, we will 
have statistically sufficient data, and the signals from GRBs will dominate the statistical 
error of ANB. Namely, then we can definitely detect the signals of the GRB neutrino back- 
ground as a significant deviation from the statistical error of ANB. The time (yr) required 
to marginally distinguish the signals from the noise is computed by t = N/S 2 , where N is 
the noise (yr -1 ) from ANB and S is the signals from GRBs (yr -1 ). In Table [III], we show 
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the time required to detect the signal of GRBs at SK/TITAND for the energy range E = 
(50 - 60) MeV and E > 50 MeV. From the table, we find that we can detect the signals 
within ~10 years at TITAND in the case that the mass-accretion rate is (5 - 10M Q s _1 ) and 
/ is (0.5 - 1). Although it seems to take a long time to detect the signals from GRBs for 
the cases in which the conservative parameter sets are adopted, constraints for optimistic 
parameter sets will be obtained by the observations of TITAND. Of course, a larger water 
Cherenkov detector than TITAND will detect the signals from GRBs more immediately even 
if the GRBs obey the conservative parameters. 

In this study, we have mainly assumed that one tenth of collapsars whose mass range 
are in (35 - 125)M are accompanied with GRBs. Here we investigate the validity of the 
assumption. The locally observed value for the SN rate is in the range -Rsn = (1-4 ± 0.5) x 
i(T 4/i o.7 yr" 1 Mpc~ 3 , where h = H / (100 km s" 1 Mpc" 1 ) f24|. On the other hand, it is 
reported that the GRB rate at z = is 2.4 x 10~ 9 /zg T yr" 1 Mpc~ 3 using the BATSE data 
catalog P5| . On the other hand, we assumed that the ratio of the collapsar rate to SN 
rate is ~ 0.13 using Eqs. (f|) and @. Thus, we estimate that the collapsar rate at z = 
is 1.8 x 10~ 5 /ip 7 yr~ x Mpc~ 3 . Adopting the beaming factor is dfl ~ An x 10~ 3 , and the 
probability that one collapsar generates a GRB is / = 0.1, we estimate the GRB rate at 
z = 0: 




R GRB = 1.8 x 10~ 9 -f- — — , (33) 



in units of yr _1 Mpc~ 3 which is consistent with the previous work. At present, the average 
beaming angle is not determined precisely except for some observations that suggest it is 
several degree pj. Therefore, we leave the value of / as a free parameter at present. In 
addition, the previous works on the GRB formation rate would have a lot of uncertainties 
since only the BATSE data catalog is used in these analyses, and the average luminosity of 
GRBs is adopted to fit the model. More precise data on the beaming angle and on the GRB 



formation rate (for example, by using the peak-lag relation |26fl) will give us informations on 
the probability /. Then, our predictions of GRB neutrino background become more correct. 

In this study, we have assumed that GRB formation rate is proportional to the star 
formation rate in the universe. If this assumption is not correct, of course, the resulting 
energy spectrum will also be changed. Therefore, on the contrary, when we observe the 
differential flux of GRB neutrino background, we will also be able to know whether the 
GRB formation history traces the star formation history or not. 

We used the model presented by Popham et al. (1999) as an accretion model. We have 
to investigate the model dependence of the neutrino luminosity in future. In particular, we 
have to investigate the neutrino-dominated accretion disk (NDAF) for the high accretion 
rate cases [0,?]- 

In this study, we have considered the effects of vacuum neutrino oscillations. However, 
the effects of MSW effects are not taken into account. As long as the collapsar model is 
adopted, the density profile of the progenitor should depend on the zenith angle and time. 
We are planning to investigate these effects in the near future. 

When the mechanism of the central engine is not the neutrino heating but the Blandford- 
Znajek mechanism ||28||, the observations of GRB neutrino background will give a lower limit 



of the total explosion energy of a GRB. 
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V. SUMMARY AND CONCLUSION 



We have estimated the flux of the GRB neutrino background and computed the event 
rate at SK/TITAND in the collapsar models, assuming that GRB formation history traces 
the star formation history. We have found that the flux and the event rate depend sensitively 
on the mass-accretion rate although the detection of signals from GRBs seems to be difficult 
by SK. On the other hand, however, we have found that the GRB neutrino background may 
be detected by TITAND within ~ 10 yrs as long as the average mass- accretion rate is high 
(> a few M s" 1 ), and the probability that one collapsar generates a GRB is high (/ = 0.5 - 
1.0). Thus, we conclude that the informations on the mass-accretion rate of collapsar will be 
obtained by the observations of the GRB neutrino background, which in turn, should give 
us informations on the total explosion energy of GRBs in future. They are, of course, never 
obtained only by the gamma-ray observations. Although there are some simplifications in 
our analyses and some uncertainties on the present observations, our proposal to determine 
the total explosion energy of a GRB is very challenging and interesting. This paper is the 
maiden attempt about the GRB neutrino background. We will revise our model to enhance 
the predictability, and our suggestions will be tested by the future observations. 
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FIG. 1. The cosmic GRB formation history. The solid, dotted and dashed curves show the 
rate of GRB formation per unit comoving volume as function of redshift for the models GRBR1, 
GRBR2, and GRBR3, respectively. A-dominated cosmology (fifcf=0.3, Q\=0.7) is adopted. The 
Hubble constant is taken to be 70 km s _1 Mpc . The probability that one collapsar generates a 
GRB (/) is set to be 0.1. 
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FIG. 2. Emissivity of V e . Here we adopt T = 5 MeV and p = 10 10 g • cm~ 3 . Solid line 
epresents the total emissivity of V e . Dashed (dotted) line is the contribution from n + e + — > p + I7 e 
e + + e _ — > f e + V e ). 
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FIG. 3. Energy spectrums of 17 e from a collapsar. Solid line represents the total energy 
spectrum. Dashed (dotted) line is the contribution from n + e + — ► p + V e (e + + e~ — > v e + F e ). 
The total accreting mass, the initial mass, and the mass-accretion rate are set to be 3OM , 3M , 
and O.IMq s , respectively. 
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FIG. 4. Duration of the neutrino emission from a collapsar as function of the mass-accretion 
rate. The total accreting mass and the initial mass are set to be 30M Q , 3M Q , respectively. 
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FIG. 5. Dependence of the energy spectrum of anti-electron neutrino on the mass-accretion 
rate (M). Solid line, dotted line, short-dashed line, long-dashed line, and dot-short-dashed line 
correspond to M = 1OM s -1 , 5M Q s _1 , 1M s' 1 , O.1M s _1 and O.O1M s _1 , respectively. 
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FIG. 6. Average neutrino luminosity as a function of M. We also show the result of Popham 
et al. (1999) for comparison. 
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FIG. 7. Efficiency of v e v e pair annihilations as a function of the mass-accretion rate, 
comparison, we also show the result of Popham et al. (1999) as a shaded region. 
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FIG. 8. Total energy of emitted neutrinos (Qtot) and total deposited energy (Q e+e _) in units 
of erg as a function of the mass- accretion rate. 
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FIG. 9. Differential number flux of GRB neutrino background [cm -2 s _1 MeV -1 ]. Solid line, 
dotted line, short-dashed line, and long-dashed line correspond to M = 10M Q s _1 , 5M s _1 , 1M Q 
s _1 , and O.IMq s _1 , respectively. The case of SN neutrino background and atmospheric neutrino 
background (ANB) are shown in dot-short-dashed line and long-short-dashed line for comparison. 
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FIG. 10. Same as Fig.9 but GRBR2 is used. 
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FIG. 11. Same as Fig.9 but GRBR3 is used. 
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FIG. 12. Event rate of GRB neutrino background [yr _1 MeV -1 ] at Super-Kamiokande. Solid 
line, dotted line, short-dashed line, and long-dashed line correspond to M = 10M Q s _1 , 5M Q s _1 , 
IMq s _1 , and O.IMq s _1 , respectively. Event rate of SN neutrino background and atmospheric 
neutrino background (ANB) are shown in dot-short-dashed line and long-short-dashed line for 
comparison. 
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FIG. 13. Dependence of the event rates on the probability / that one collapsar generates a 
GRB. Solid line and dotted line correspond to the case / = 1.0 and / = 0.1, respectively. We 
adopt the case of M = IOMqs" 1 
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FIG. 14. Event rate at SK for the model of M = IOMq. Solid lines represent the case in 
which the effects of neutrino oscillation are taken into consideration. Dotted lines represent the 
case of no oscillation. 
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TABLES 



TABLE I. Event Rate (yr^ 1 ) and Statistical Error of ANB at SK/TITAND 

N (yr- 1 ) VNU^) 
E Pe = (50 - 60)MeV 1.0/1.0xl0 2 1.0/l.OxlO 1 

E Pe > 50MeV 2.6/2.6xl0 2 1.6/1.6X10 1 



TABLE II. Event Rate (yr^ 1 ) of Signals at SK/TITAND 



(Mqs- 1 ) 


/ 


Oscillation 


GRBR 


E Pe =(50 - 60)MeV 


E Pe >50MeV 


10 


0.1 


Yes 


3 


9.1xl0~ 3 /9.1xl0 _1 


1.7xl0" 2 /1.7 


10 


0.5 


Yes 


3 


4.5xl0~ 2 /4.5 


8.6xl0- 2 /8.6 


10 


1.0 


Yes 


3 


9.1xl0- 2 /9.1 
1.8x10 z /1.8 
9.0xl0- 2 /9.0 


1. 7x10-71. 7X10 1 
3.4x10 Z /3A 


10 


0.1 


No 


3 


10 


0.5 


No 


3 


1.7xl0-7 1 - 7xl ° 1 


10 


1.0 


No 


3 


1.8X10-7L8X10 1 


3.4x10-73. 4X10 1 


5 


0.1 


Yes 


3 


4.6xl0" 3 /4.6xl0" 1 
2.2xl0" 2 /2.2 
4.6x 10~ 2 /4.6 


8.3x10-78. 3X10- 1 


5 


0.5 


Yes 


3 


4.2xl0" 2 /4.2 
8.3xl0- 2 /8.3 
1.7xl0- 2 /1.7 


5 


1.0 


Yes 


3 


5 


0.1 


No 


3 


9. lxl0- 3 /9. lxlO- 1 


5 


0.5 


No 


3 


4.6xl0- 2 /4.6 


8.3xl0- 2 /8.3 


5 


1.0 


No 


3 


9.1xl0- 2 /9.2 


1.7xl0 _ 7 1 - 7xl ° 1 




0.1 


Yes 


3 


8.3xl0-78.3xl0- 2 


1.4xl0- 3 /1.4xl0 _1 




0.5 


Yes 


3 


4.2X10-74.2X10- 1 


6.9X10-76.9X10- 1 
1.4xl0 -2 /1.4 




1.0 


Yes 


3 


8. 3x10-78. 3X10- 1 




0.1 


No 


3 


1.7xl0 -3 /1.7xl0 _1 


2. 7x10-72. 7X10- 1 




0.5 


No 


3 


8. 3x10-78. 3X10- 1 
1.7xl0 -2 /1.7 


1.4xl0 _2 /1.4 
2.7xl0~ 2 /2.7 




1.0 


No 


3 
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TABLE III. Time (yr) Required to Detect the Signals at SK/TITAND 



v \y / 


/ 


Oscillation 


GRBR 


E 9 =(50 - 60)MeV 


E p >50MeV 


10 


0.1 


Yes 


3 


1.3xl0 4 /1.3xl0 2 


2.3xl0 4 /2.3xl0 2 


10 


0.5 


Yes 


3 


5.1xl0 2 /5.1 


9.3xl0 2 /9.3 


10 


1.0 


Yes 


3 


1.3xl0 2 /1.3 


2.3xl0 2 /2.3 


10 


0.1 


No 


3 


s.ixioVs.ixio 1 

1.3xl0 2 /1.3 


5.8X103/5.8X10 1 
2.3xl0 2 /2.3 


10 


0.5 


No 


3 


10 


1.0 


No 


3 


3.2X10V3.2X10" 1 


5.8xlOV5.8xlO- 


5 


0.1 


Yes 


3 


5.0xl0 4 /5.0xl0 2 


9.9xl0 4 /9.9xl0 2 


5 


0.5 


Yes 


3 


2.0xl0 3 /2.0xl0 1 
5.0xl0 2 /5.0 


4.0xl0 3 /4.0xl0 1 
9.9xl0 2 /9.9 


5 


1.0 


Yes 


3 


5 


0.1 


No 


3 


1.3xl0 4 /1.3xl0 2 
5.0xl0 2 /5.0 


2.5xl0 4 /2.5xl0 2 

9.9xl0 2 /9.9 

2.5xl0 2 /2.5 


5 


0.5 


No 


3 


5 


1.0 


No 


3 


1.3xl0 2 /1.3 


1 


0.1 


Yes 


3 


1.5xl0 6 /L5xl0 4 


3.7xl0 6 /3.7xl0 4 


1 


0.5 


Yes 


3 


6.0xl0 4 /6.0xl0 2 


1.5xl0 5 /L5xl0 3 


1 


1.0 


Yes 


3 


1.5xl0 4 /1.5xl0 2 


3.7xl0 4 /3.7xl0 2 


1 


0.1 


No 


3 


3.8xl0 5 /3.8xl0 3 


9.2xl0 5 /9.2xl0 3 


1 


0.5 


No 


3 


1.5xl0 4 /1.5xl0 2 


3.7xl0 4 /3.7xl0 2 


1 


1.0 


No 


3 


s.sxio-Vs.sxio 1 


9.2xl0 3 /9.1xl0 1 
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